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We study the effect of including a BGK coUisional kernel on the collective modes of a QCD 
plasma which has a hard-particle distribution function which is anisotropic in momentum space. We 
calculate dispersion relations for both the stable and unstable modes and show that the addition of 
hard particle collisions slows the rate of growth of QCD plasma unstable modes. We also show that 
JH ■ for any anisotropy there is an upper limit on the coUisional frequency beyond which no instabilities 

' exist. Estimating a realistic value for the coUisional frequency for ~ 0.2 — 0.4 we find that for the 

large-anisotropy case which is relevant for the initial state of matter generated by free streaming in 
heavy-ion collisions that the coUisional frequency is below this critical value. 
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I. INTRODUCTION AND MOTIVATION 
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The ultrarelativistic heavy ion collision experiments ongoing at the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory (BNL) and planned at the Large Hadron Collider (LHC) at CERN study nuclear 
matter under extreme conditions. A main point of interest is the identification and investigation of a phase transition 
to the quark-gluon plasma (QGP). If it is created, the QGP is expected to expand, cool and finally hadronize. Whether 
I or not it is possible to apply a thermodynamic description of the created system in its later stages of evolution depends 
^p^i on whether it thermalizes fast enough. RHIC data indicate that such a rapid thermalization does occur 0,0,13' 
O 1' contradiction to leading order perturbation theory estimates. Therefore, many recent works have been dedicated to the 
, explanation of how the fast thermalization is achieved. One possibility is the assumption of a strongly coupled QGP 
' • In other approaches @ the process of particle production leads to momentum distributions of the equilibrium 

^ form immediately, without any secondary processes needed. However it is not explained how the equilibrium state is 
maintained when the system is driven out of equilibrium by free streaming. Secondary processes are certainly necessary 
to explain this. Within recent transport theory approaches, the inclusion of particle production and absorption via 
2 <-> 3 pQCD bremsstrahlung processes speeds up the equilibration significantly as compared to equilibration solely 
driven by perturbative binary collisions. 

Here we will concentrate on the role of non-equilibriurn QCD collective modes in the isotropization and thermal- 
ization of the system as investigated in |1 11 [Tfl, [ll [ll [ll [ll [ll [ll [ll [ll [ll |M IM (for a recent review 
see H^l). In these approaches the QGP is assumed to be homogeneous and stationary, but anisotropic in momentum 
space, such that kinetic instabilities can occur. They are initiated either by charge or current fluctuations. In the 
first case, the electric field is parallel to the wave vector k (E || k), while in the second case the field is perpendicular 
to k (E _L k). This is why the corresponding instabilities are called longitudinal and transverse, respectively. Since 
the electric field plays a crucial role in the generation of longitudinal modes, they are also called electric, while the 
transverse modes are called magnetic. The magnetic mode known as the filamentation- or Weibel-instability 
appears to be relevant for the QGP [1,0,01, created in heavy ion collisions. 

The instabilities are found to have a significant effect on the system's evolution, leading to a faster isotropiza- 
tion and equilibration. The equilibration due to instabilities only happens indirectly, because the instabilities driven 
isotropization is a mean-field reversible process, which does not produce entropy. However, parton momentum dis- 
tributions are influenced by the isotropization, which speeds up the equilibration. Collisions, being responsible for 
the dissipation are needed to reach the equilibrium state of maximum entropy. All calculations so far have been 
carried out at leading order in perturbation theory, such that collisions among the hard particles that enter at higher 
orders in g could be neglected. However, in heavy ion collision experiments at the RHIC and LHC the couplings 
expected are of the order as ^ 0.2 — 0.4 and higher order terms will be important. Hence, collisions can not simply 
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be neglected and their effect on the system's evolution, particularly on the collective modes, has to be investigated. 
We expect that the inclusion of coUisional damping will decrease the instability growth rate, however, one would like 
to make this statement quantitative and say precisely by how much the growth rate is reduced. This is the main 
objective of the present work in that we give a first quantitative estimate of how the collisions affect the dispersion 
relations of the collective modes and the growth of instabilities in particular. To achieve this, we introduce a model 
for the inclusion of collisions based on the Vlasov equations for QCD combined with a BGK-type collision term. 
The regarded anisotropic distributions are derived from an isotropic one by contracting the hard particle distribution 
function along a preferred direction. We concentrate on the case in which the direction of anisotropy is parallel to 
the wave vector of the regarded collective mode, because in this case the growth rate of the magnetic instability is 
maximal and analytic expressions for the structure functions can be found. For any such distribution we find that the 
growth rate of instabilities decreases approximately linearly with increasing collision rate and there exists a critical 
collision rate above which instabilities cease to exist. 

The current work is organized as follows: In Section ^] we briefly review the kinetic approach in the coUisionless 
case. We introduce the model for the inclusion of collisions in Section UTTI which closes with the final result for the 
self energy of the collective modes. Section Hvl bricflv reviews the procedure for finding the dispersion relations and 
in Section we derive the structure functions and dispersion relations for the specified anisotropic systems. Final 
results for the dispersion relations of the stable modes including collisions are given in Section [VII while the influence 
of collisions on the unstable modes is shown in Section lVlII In Section lVlIII we discuss further the results and attempt 
to estimate the numerical value of the coUisional frequency. In Section Hxl we conclude and give an outlook. 



II. COLLISION FREE TRANSPORT EQUATIONS 



In the kinetic theory approach |^ 0, 0, HE HE HE HEl gluons, quarks and antiquarks are described by 
their phase space densities, given by the gauge covariant Wigner functions X), with i e {g, q, q}. To obtain the 
linearized transport equations we expand n^{p,X) = n^{p) + Sn''{p,X). After a gauge covariant gradient expansion 
of the generalized QCD-Kadanoff-Baym equations for the quark and gluon Green functions, the linearized equations 
of motion read |1 i, [IE 111 111 HI ll HI 111 113] 



[V ■ Dx,Sn\p, X)] + g0,V^F^''{X)djj''>n\p) - , 



(1) 



where V — (1, v) with v = p/|p| and dg = 9q = 1 and 9q = —1. Dx ~ dx + igA{X) is the covariant derivative with 
the soft gauge flcld = A^T" or = A^^t"- in the gluon and quark equations of motion, respectively. The gluon 
field strength tensor is given by F^'^ — d^A'^ — d'^A^^ — ig[A^^ , A'^]. Eq. ^ holds for each particle species, with the 
color neutral background fields n'^^^ij)) — /'^'(p)/ and n^{p) = /^(p)X. / andX are unit matrices in the fundamental 
and adjoint representation, respectively. The scalar functions (p) and /^(p) arc found by projection: 
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r/Hp) = —TT[n'^/Hp,X) 



m-1 



TrK(p,X)] 



(2) 



The induced fiuctuations = 6f^'^{p,X)t^ and 5n3{p^X) = 6fjj{p,X)T\ with the generators U and 

Ti, in the fundamental and adjoint representation, respectively, are assumed to be much smaller than the colorless 
background terms. The scalar functions dfl{p,X) again follow by projection: 



Sf^/\p,X) = 2TT\tan^/\p,X) 



Sf^{p,X)^—TT[Tan3(p,X)] . 



(3) 



The induced current for each color channel reads [ill H3| 



JZa ^9 {wmp, X) + Ni[5n{p, X) - sriip, x)]} , 



(4) 



with 



(2^ 



(5) 
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When we concentrate on the soft scale k ~ gT <C T, the first scale at which collective motion appears, we can neglect 
terms of subleading order in g and the theory becomes effectively Abelian.^ Now all color channels decouple in Eq. 
and we find separate solutions for each Sfl{p, X). After Fourier transformation the total induced current Q becomes 

m 

JLaiK) - g' V%^fip) (g,, - /'yl^^ ) AJiK) + 0{g'A') , (6) 

where K =^ ("-^i k) and the effective background phase space density is given by 

/(p) = 2A^,/^(p) + Nf [Pip) + fip)] . (7) 

Functional differentiation of the induced current with respect to the gauge field and taking ^ leads to the 
gluon self-energy, which equals the one which can be obtained by a hard-loop resummation within the diagrammatic 
approach in the coUisionless case. 



III. INCLUSION OF COLLISIONS 

In the effectively Abelian limit, introduced in Section^ the equation of motion holds for each color channel 
separately, such that we can write 

V ■ dxSflip, X) + ge.V^Fi:%X)d^J^f\p) = Clip, X) , (8) 

where we have included the collision term C^, given by 



Cl{p,X)^~iy 



fl{p,X)~^^f,^{\p\) 



(9) 



with fl{p,X) = /'(p) + 6fl{p, X). This BGK-type collision term [23 describes how collisions equilibrate the system 
within a time proportional to h'~^. Here we will assume that the collision rate u is independent of momentum and 
particle species, however, these assumptions are easily relaxed. Note that these collisions are not color-rotating. The 
particle numbers are given by 

Kix)= f ra{p,x) , f /<!q(ipi)- / r(p). (10) 

We note that the difference between the BGK coUisional kernel @ and the conventional relaxation-time approximation 
(RTA) is the multiplication of the second term in brackets by the ratio of the density over the equilibrium density. RTA 
simply takes the difference of the distribution function and the equilibrium distribution function implicitly setting 
this ratio to one. The advantage of the BGK kernel over an RTA kernel lies in the fact that the number of particles 
is instantaneously conserved by the BGK coUisional- kernel, i.e., that 

c:(p,x) = o. (11) 



p 



This simply states that the collisions can only occur if a particle is present and that only the momentum of the 
particles will change as a result, not the particle number. This condition is violated by RTA. In addition, Ref. [sJl 
showed that after a simple modification the BGK coUisional- kernel covariantly conserves the color current. 

The inclusion of a BGK coUisional kernel allows for simulation of the effect of binary collisions with substantial 
momentum transfer and is merely an approximation for collisions between the hard charged particles in a hot quark- 
gluon plasma. However, in our opinion it is a reasonable way to yield a first quantitative answer to the question of 
how collisions among hard particles affect the collective modes of QCD. The effects of such a collision term on the 
dispersion relations in the ultrarelativistic case for an isotropic system were investigated in [s^ . For now the collision 



^ This strictly only occurs at leading order in g. The induced current also contains terms of higher order in g (and A), which correspond 
to the non- Abelian self-interactions of the soft gauge field. Here we are looking at the self-energy, so that we can ignore these. 



rate v is taken to be a free parameter and we postpone the estimation of its magnitude to the discussions at the end 
of this paper. 

Using we can write 



r(p) + <5/:(p,x)- 1 + 



/eq(IPl) 



cq 



(12) 



Solving for 5p^{p, X) and Fourier-transforming leads to the result for the linearized induced current by each particle 
species i (see Appendix El : 



with 



+ 9^ / V^/e\(|p|)l?-^(if,V,i.) 
cq -/ p 



X 9 



p' 



dl,^r{p')M,p{K, V')D-\K, v', i^)Al{K) + giyS\K, v) 



Mjp{K, V) := g^p{uj - k • v) - V^Kfs , 



(13) 



(14) 



D{K, V, i^) := u! + iu — k ■ V , 



(15) 



and 



S\K,,.) -.^9, f y^[r(p)-/^q(|p|)p-i(if,v,, 

Jp 

W:{K,v) :=1--^ / f:^{\Yy\)D-\K,w,v). 



(16) 



(17) 



The total induced current is given by Jf^^^jK) = 2NcJ[-l^^iK) + Nf [J^^'^JK) + J^^'^JK)] . It can be simplified due 
to the fact that the integrals over f^^ in Eq. 1131) are independent of the particle species: 



(18) 



where dfl = sin 9d9dip. Assuming equal distributions for quarks and antiquarks, the full linearized induced current 
reads 



JLaiK)=9' I V'^d^^^f{p)M^i,{K,V)D-\K,v,,.)A2+2N,g,.S^{K,u) 



p' 



dVL^ 



+ 2Ncg^(iiy^) / — V''D-^(K,yr,iy)S3(K,iy)W-\K,v) 



(19) 



where W{K, ly) — 1 — iv J ^{K, v, v) and /(p) as in Eq. 0. The self energy is obtained from Eq. (fT^ via 



(20) 
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resulting in 

n^rW^ W / V^df f{p)M''^{K,V)D-\K,v,i.) 
Jp 

+ 5ab9\iy) I ^V>^D~\K,v,i.) j d'-p f{v')M''^{KX)D~\Ky ,v)W-'{K,v) , (21) 

which is diagonal in color and can be shown to be transverse, i.e., K^H'^'^ = Ki^H^" = 0. We will from now on 
omit the color indices of 11^'' . The terms in the induced current involving 5^ drive the distribution into an isotropic 
equilibrium shape. They represent a parity conserving current and thus create a zero average electromagnetic field 
and therefore does not contribute to instability growth. 



IV. DISPERSION RELATIONS 



In the linear approximation, the current that is induced by the fluctuations can be expressed in terms of the self 
energy by 

JL^{K)^n^''{K)A,{K), (22) 

where Jf^^ = Jf^^ ^ + Jl^-^^ 2 ^^'^ "^hid 2 ^^'^ parity conserving part of the current l(T^ that does not couple to the 
gauge field. Note that for a parity symmetric distribution as we use below, J/^^jj ~ ^- Inserting this into Maxwell's 
equation 

iif^F^''(X) = Ji"^di + ^cxt, (23) 

we obtain 

[K^g^'' -K^K'' + n^"^ {K)]A^{K) = J^^, [K) , (24) 
with the external current J^xt- In the temporal axial gauge, where = 0, this becomes 

[(fc2 - _ 4, IT^K)] E^K) = [A-\K)f (K) ii^r^^^K) , (25) 

and the response of the system to the external source is given by 

E^iK)^ZLuA^^iK)JUK). (26) 
The dispersion relations are obtained by finding the poles of the propagator A^^{K). 



V. SELF ENERGY AND PROPAGATOR IN AN ANISOTROPIC SYSTEM 

So far, the distribution function /(p) in Eq. (|21|) has not yet been specified. We assume that /(p ) ca n be obtained 
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from any isotropic distribution function by rescaling one direction in momentum space by defining 

/(P) = yr+^/iso (p' + e(p • fi)) , (27) 

for an arbitrary isotropic distribution function /iso(|p|)- Distributions like (|27|l with ^ > are generated during 
a heavy-ion collision at times t > {pt)~^. The direction of the anisotropy above is given by n and ^ > — 1 is an 
adjustable anisotropy parameter. ^ > corresponds to a contraction of the distribution in the n direction, whereas 
— 1 < ^ < represents a stretching of the distribution in the n direction. The factor of vT+T ensures that the overall 
particle number is the same for both the anisotropic and the isotropic distribution function. With this particular 
distribution we are able to perform the radial part of the integrations involving /(p) in Eq. (|21|) by changing variables 
to = p^{l + ^(v • n)2). The result is 
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where 



dp 



(29) 



We now decompose the self energy into four structure functions, using the general tensor basis for an anisotropic 
system, developed in 0|, such that 



where 



(30) 



(31) 



with — A^^n^ the part of n that is perpendicular to fc, i.e., fi ■ k — Q. We determine the structure functions by 
taking the contractions: 



k'W^y = n'ffJp = fi^k^S, nm'^n^ = ^^(a + 7), Trff^' = 2a + /3 + 7 . 
In this basis, the inverse of the propagator in Eq. H25I) can be written as 

A-\K) = (fc2 -uj^ + a)A + (/3 - io^)B + jC + SB, 
while the propagator itself is given by p^ : 

A(/^) = Aa{A - C) + Ag [(fc' - + a + 7)B + (/? - cj2)C - 5D] , 

with 

A^^{K) = k^-cj^+a, 

A^^{K) = (fc2 - + a + 7)(/3 - w^) „ ^2-2^2 ^ 



(32) 

(33) 

(34) 

(35) 
(36) 



Since the growth rate of the filamentation instability is the largest when the wave vector is parallel to the direction 
of the anisotropy, i.e., k || n (see e.g. ^3,^Sl23|)i '^^ concentrate on this particular case. Then 



k-v = fcn-v = fc< 



(37) 



and 7 and — 1 — (k- n)^ vanish identically. To determine the poles of the propagator (|34|l . and hence the dispersion 
relations, we are now left with two separate equations for the a— and /3— mode, respectively: 



yfc2 - _^ a = , 
/3 - = . 



(38) 



From Eqs. (|32|l we find a and /3, which correspond to Ht and (w^/fc^)!!^, respectively. The integrals can be solved 
analytically and the final results simplify to: 



a{uj,k,^,u) 



4 A:(l+ez2)2 
i 



{k{z^ - 1) - iziy)il + ^z^) - (z^ - l)(fcz(l + - it") In 



z + 1 



z-1 



zi^{l + {3 + z\l- 0)0 + - e + 2^(1 + e(6 + z2(e - 1) + 0))] arctan ^1 \ , 



Pico, k, ^) = nil.JT+lk{kz - ivf I - 271(1 + 7?0 + (1 + (2z 



■In 



z + 1 
z - 1 



+ 2(z2$- 1) arctan 



2v/^(l + z^O^fc^ (2fc-ii/ln 
where we abbreviate z = (cj + w)/k. 



z + 1 
z-l 



(39) 



(40) 
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FIG. 1; Real part of the dispersion relation for the stable Q-mode for an anisotropy parameter of ^ = {0, 1, 10} and different 
collision rates in units of mg = mn / VS. 
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FIG. 2: Imaginary part of the dispersion relation for the stable a-mode for an anisotropy parameter of ^ = {0, 1, 10} and 
different collision rates in units of rug = mo/VS. 



VI. STABLE MODES 



The dispersion relations for all modes are given by the solutions uj{k) of Eqs. H38(l . These solutions are found 
numerically for different values of the collision rate v. The results for the stable transverse (a-) mode are shown 
in Figs, n and |21 for two anisotropic distributions, ^ = 1 and ^ = 10, together with the result for the isotropic case 

= 0). Note that the finite collision rate causes uja to become complex with a negative imaginary part corresponding 
to damping of these types of modes. 

The effect of collisions on the stable longitudinal (/?-) modes is more significant. The results are presented in Figs. El 
and^l We find that for finite v the dispersion becomes spacelike (Re(w) < k) at large k in contrast to the coUisionless 
case, in which Re(ti;) > k always holds. This behavior is responsible for the fact that the solution vanishes from the 
physical Riemann sheet above some finite k. This occurs precisely when the solution for u!/j reaches the cut between 
—k and k at —iv. 

It is however possible to find solutions on the unphysical Riemann sheets by replacing the logarithm in the structure 
function with its usual analytic continuation [T^ : 



In 



z + 1 



= In 



z + 1 



+ i 



arg 



z + 1 



2nN 



(41) 



where N specifies the sheet number. The continuation of the solution to the lower Riemann sheets is shown in Figs.^l 
and for a collision rate of v = 0.8 « 0.46 mu and different anisotropy parameters ^. For smaller anisotropics 
the solution is found to converge to the light cone in an oscillating manner, while the imaginary part of lu oscillates 
around +v. Between ^ = 2 and C = 3 (for v w 0.46 m^i) this behavior changes qualitatively to the one shown in the 
case ^ = 10. With increasing k the real part of u moves away from the light cone, while the imaginary part drops to 
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FIG. 3; Real part of the dispersion relation for the stable /3-mode for an anisotropy parameter of ^ = {0, 1, 10} and different 
collision rates in units of rUg = mo/VS. Note that for finite v there is a maximal k beyond which there is no solution. It 
vanishes when uj crosses the cut in the complex plane, which extends from ~k to +k at ~iu (cf. Fig|l|l. 
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FIG. 4: Imaginary part of the dispersion relation for the stable /3-mode for an anisotropy parameter of ^ = {0, 1, 10} and 
different collision rates in units of nig = mo/VS- The solution vanishes, when uj crosses the cut in the complex plane, which 
extends from — fc to +k at —iv (indicated by the straight dotted line). 



large negative values with the solution remaining on the = — 1 Riemann sheet, such that these modes are strongly 
damped. 



VII. UNSTABLE MODES 



Anisotropic momentum distributions cause kinetic instabilities. For the particular distribution (j27|l with n || k 
and ^ > there exists a magnetic instability, the so called filamentation or Weibel instability [l^. Its existence is 
due to a surplus of particles with momentum perpendicular (or close to perpendicular) to k. These particles are 
trapped in the direction of k by the background magnetic field and cause currents, which generate magnetic fields 
that add to the original one. Hence they contribute to instability, while all other particles have a stabilizing effect, 
(for a more detailed qualitative review of this scenario see [Tsl ) . In the isotropic case the stabilizing and destabilizing 
contributions cancel, such that no instability arises. 

We now investigate how the inclusion of collisions as described in Section IIIII affects the growth rates of these 
instabilities. Qualitatively one expects a decrease of the growth rates because the particles, which move perpendicular 
to k, can scatter with other particles and will no longer be trapped. Other particles can gain a momentum close 
to perpendicular to k and form a new contribution to the instability. However, since the collision term tends to 
randomize the momentum distribution, the growth of Sf and the magnetic field is prevented. In order to describe this 
effect quantitatively, we solve Eq. I|38|) for purely imaginary uj and vary the collision rate v. The solution uj{k) = iT{k) 
gives the growth rate r(fc). In the case that n || k solutions like that only exist for the transverse (a-) mode. The 
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FIG. 5: Real part of the dispersion relation for the stable /3-mode for an anisotropy parameter of ^ = {0, 1, 10} and u — 0.8 mg. 
It is shown how the solution continues on lower Riemann sheets. 
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FIG. 6: Imaginary part of the dispersion relation for the stable /3-mode for an anisotropy parameter of ^ = {0, 1, 10} and 
v — 0.8 mg. It is shown how the solution continues on lower Riemann sheets. Note the different scale for the third plot. 



one with F > corresponds to the filamcntation instability. Results for different values of the collision rate v are 
shown in Fig. |7|for ^ — I and in Fig. [S] for ^ = 10. The qualitatively expected effect is nicely reproduced. The growth 
rate decreases with an increasing collision rate as does the maximal wave number for an unstable mode. One can 
see that in the case where ^ = 1 already for i/ being around 20% of the Debye mass, growth has completely turned 
into damping and no instability can evolve. For ^ = 10 the collision rate v has to be slightly larger than 30% of the 
Debye mass in order to prevent growth of a collective mode. In order to find the wave number A;max(^, i^) at which 
the unstable mode spectrum terminates, we take the limit lu ^ to obtain 



ml 



-2t/(fc2 + t.2)^3/2 _ -^r^) - 2«fc - 1) + + 3)] arctan(v/C) | 



(42) 



One of the solutions to the equation 



k'+mi = 0, (43) 

which is just the limit ^ of the first of the Eqs. H38|l . is fcmax- Results for different ^ and i/ are shown in Figs.|51 
and 1 101 We find that for a given anisotropy parameter ^, there exists a critical collision rate, above which instabilities 
can not occur. This is also true for the limit ^ oo, as we will show in the following. Taking ^ oo, Eq. 1)42(1 
becomes 

TT P 

mli^^oo) = --ml-^, (44) 
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FIG. 7: Dependence of the growth rate F of the unsta- 
ble transverse (a—) mode on the colhsion rate i/, for an 
anisotropy parameter ^ = 1. 



FIG. 8: Dependence of the growth rate F of the unsta- 
ble transverse (a—) mode on the collision rate i/, for an 
anisotropy parameter ^ = 10. 
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FIG. 9: fcniax of the maximally unstable mode as a func- 
tion of the anisotropy parameter ^ for different values of v. 
For given v, the momentum anisotropy parameter must 
be above the critical value for instabilities to occur. 



FIG. 10: fcmax of the maximally unstable mode as a func- 
tion of the collision rate v for given anisotropy parameter 
^. There exist ^-dependent critical collision rates !-'max(C) 
above which instabilities can not exist. 



which together with Eq. H43|l gives 

k\,.'-^rnl)=Q. (45) 

Apart from the result fe = 0, this is solved by 

t'maxCC —^oo) = ^J^rno « 0.6267 TOD = T^axlC °°) - (46) 

the critical collision rate, above which even in the extremely anisotropic limit ^ = oo no instability can occur. 

This is also visible in the plot of the growth rate T^^oo, shown in Fig. 1111 For v = z/max(C oo) the growth rate 
becomes zero for all k, and for larger v only damping occurs. In this case = oo) the value of the maximal collision 
rate equals that of the maximal growth rate in the coUisionless limit. This simply means that the instability vanishes 
completely at the point where the collisions damp at the same rate at which the instability grows. Note however that 
this relation is more complicated in general as shown in Fig. 1121 where the dependence of the maximal growth rate on 
the collision rate is plotted. In order to make the instability vanish completely for any ^ < oo, a. collision rate larger 
than the maximal growth rate of the instability in the coUisionless limit is needed. 



VIII. DISCUSSIONS 



In the body of the text the coUisional frequency v has been taken to be arbitrary. Because the inclusion of a 
BGK coUisional term is a phenomenological model for the equilibration of a system and cannot be derived from 
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FIG. 11: Dependence of the growth rate F of the unstable transverse (a—) mode on the colhsion rate v, for the extremely 
anisotropic limit ^ = oo. 
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FIG. 12: The maximal growth rate of the instability as a function of the collision rate u. 



first principles this makes it difRcult to fix the magnitude of u. However, since as treated here the underlying 
framework of Boltzmann-Vlasov is implicitly perturbative, we can attempt to fix v perturbatively. However, even this 
is non-trivial since within non-abelian theories there are at least two possible collisional frequencies to be considered 
H^ISHEUE 37J: (1) the frequency for hard-hard scatterings which is parametrically z^hard ^ a1 logaj^ and (2) the 
frequency for hard-soft scattering which is parametrically z/goft ~ Q^s logaj^. 

The hard-hard scatterings correspond to interactions which change the momentum of a hard particle by ©(phard) 
and therefore represent truly momentum-space isotropizing interactions. On the other hand the hard-soft scatterings 
correspond to changes in momentum which are order ©(gphaid)- These small angle scatterings occur rather frequently 
and it turns out that after traversing one hard scattering mean free path, Ahard ~ '^hard' ^he typical deflection of the 
particle is also 0(1).^ The physics of small- angle scattering by the soft-background is precisely what is captured by 
the hard-loop treatment. However, the hard-loop framework doesn't explicitly take into account that Vsoit is also 
the frequency at which there are color-rotating interactions of the hard particles themselves. One would expect that 
color-rotation of the hard particles to have a larger effect on the growth of instabilities than the momentum-space 
isotropization via hard-hard scattering. That being said, the form of the BGK scattering kernel does not mix color 
channels and in that sense cannot be used to describe the physics of color-rotation of the hard particles. For this reason 
one is lead to the conclusion that when using the BGK kernel the appropriate damping rate is ~ t'hard ~ oii log oi~^ . 
This conclusion for the parametric dependence of v is also supported by looking at the leading order result for the 



^ This, in the end, is the source of the logarithm in I'hard above. 
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shear viscosity [Sg . 

Even with this conclusion it is hard to say anything quantitative about v since the overall coefficient and the 
coefficient in the logarithm are not specified by such a parametric relation. One could hope that previous calculations 
contained in of parton interaction rates could be of some use. Unfortunately, for a purely gluonic plasma it was 
found V — 5.2 a^T log(0.25Q!7^), which clearly cannot be trusted for the values of which are relevant for heavy-ion 
collisions {as ~ 0.2 — 0.4) since i' becomes negative for large as- We note here that the negativity of this result at 
large couplings most likely stems from the strict perturbative expansion of the integrals involved failing when the hard 
and soft scales become comparable in magnitude. Similar erroneous negative values also occur in the perturbative 
expressions for heavy-quark coUisional energy loss when extrapolated to large coupling. A corrected calculational 
method which yields positive-definite results for the heavy-quark energy loss was detailed in Refs. |40ll4ll|. 

Ideally, one would revisit the calculation of the interaction rate and improve upon the techniques used where 
necessary. Short of such a calculation one cannot say with certainty what the numerical value of v should be and the 
best we can do is to play "games" . For instance, one could insert a one into the logarithm appearing in similar 
to what other authors have done to obtain ^ ^ 5.2asTlog(l -I- 0.25aJ^) ^ 0.1 - 0.2mD for as = 0.2 - 0.4. 

1 /2 

Using this admittedly specious expression for large-coupling one can expand and obtain limc^^oo v/tud ~ 0.37 . 
Note that in the case Nf = 2 both v and mo increase; however, the ratio of these two scales is still in the range 
quoted for Nf — 0. The range, v = 0.1 — 0.2 mu, places us well below the threshold needed to turn off instabilities 
in the case of extremely anisotropic distribution functions but does imply (see Fig. I12|l that for moderately small 
anisotropies, C ~ Ij ^^^d large coupling that it is possible for coUisional damping to eliminate the unstable modes from 
the spectrum completely. Of course, in the limit of asymptotically small couplings the ratio v/rriD approaches zero 
and the coUisionlcss results hold to very good approximation. In the opposite limit of strong coupling the estimates 
here are at best guesswork and it is inde ed p ossible that the ratio v /rriD is larger than the range we have quoted. For 
example, the recent work of Peshier [43ll44| implies that v/mo could be as large as 0.5; however, this number results 
from a fit of a model assumption to lattice data and is not directly comparable to the coUisional widths considered 
here since in their description the gluon width was assumed to be parametrically given by Vsoft- 

Additionally, we have to mention the caveat that all the estimates above rely on full equilibrium thermal field 
theory calculations. For the very initial state of the matter created in an ultrarelativistic heavy-ion collision the 
system is clearly not in equilibrium and it is not clear how this estimate will change as a result. However, it is of 
crucial importance to attempt to estimate the scattering rate in a non-equilibrium setting. 



IX. CONCLUSIONS 



In this paper we have studied the effects of including a BGK collision kernel on the collective modes of a QCD 
plasma which has a hard particle distribution function which is anisotropic in momentum space. To simplify the 
analysis we have specialized to gluonic collective modes which have their momentum vector, k, directed along the 
anisotropy direction, n. The reasons for doing this were two-fold: (1) in the coUisionless limit these modes correspond 
to the modes which arc the most unstable and (2) when k || n it is possible to obtain analytic expressions, Eqs. I|39|) 
and (|40|l . for the soft-gluon self-energy structure functions. We have presented herein dispersion relations for both 
the stable and unstable modes in the case that there is a finite coUisional frequency (or damping coefficient), v. Our 
results confirm what can be expected intuitively, namely that the addition of coUisional damping slows down the rate 
of growth of the unstable modes. 

However, going beyond this intuitive expectation we have presented detailed calculations of the dependence of 
the maximal unstable mode growth rate on the parameter v as shown in Fig. 1121 For all values of the anisotropy 
parameter, ^, we find that there is a critical value of i^max above which no instability is present in the system. In 
Fig. 1121 this corresponds to the value of i' at which the maximal growth rate vanishes. In the limit that ^ — > cxd we 
were able to derive an analytic expression for Vmax, Eq. H46() . finding that it corresponds precisely to the maximal 
growth rate obtained in the coUisionless limit. 

In addition, we have investigated the non-trivial analytic structure of the soft-gluon propagator in this model finding 
that the stable longitudinal mode becomes spacelike'^ and only remains on the physical Riemann sheet up to a certain 
critical momentum. Beyond this critical momentum the solution goes through the logarithmic cut to the iV = — 1 
Riemann sheet. For weak-damping the longitudinal mode then continues to spiral around the logarithmic branch 
point onto lower and lower Riemann sheets as its momentum increases (see Figs. and|S^ and Figs and|^). 
In the case of stronger anisotropies the longitudinal mode solution no longer "spirals down" the logarithmic branch 



^ In the sense that Re{Lij^/fc) < 1. 
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point to lower Riemann sheets but instead simply moves lower in the complex plane of the N = —I Riemann sheet 
(see Figs. andEt). We note that even in the isotropic case, the detailed analytic structure of the coUisionally 
damped modes enables one to calculate quantities such as the QCD pressure including the effects of damping of the 
quasiparticle modes which could be of some interest psl IS liM IIgL I47L |43 | . 

In the discussions section of the manuscript we have attempted to fix a numerical value for the coUisional frequency 
ly. We have argued that since the BGK kernel does not rotate the color of the hard particles that the appropriate 
frequency is parametrically given by the time scale for hard-hard collisions, namely v a^logo:^^. Going further 
than this parametric estimate to a value applicable at couplings expected to be generated during heavy-ion collisions 
{as ^ 0.2 — 0.4) is problematic due to the small coefficient which appears in the logarithm resulting in v becoming 
negative. Playing a game by adding a one in the argument of the logarithm we found that v ~ 0.1 — 0.2mjj which 
according to the results of this paper would imply that for weak anisotropics there are no instabilities (see Fig. I12|l . 
For stronger anisotropics the results contained here tell us how much the maximal growth rate of the unstable modes 
is affected by the inclusion of coUisional damping via a BGK kernel. Looking forward more detailed calculations of 
collision rates in a time-evolving soft-field background and true non-equilibrium situation are clearly needed. 
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APPENDIX A: ANALYTICAL SOLUTION OF THE LINEARIZED TRANSPORT EQUATIONS 

We briefly sketch how the result Eq. H13|l emerges from the transport equation (|12|l . From Eq. (|12|l we immediately 
get 



/cq(P) 

^cq Jp' 

where 6f^{p,K) and F'^''{K) are the Fourier-transforms of Sf^{p,X) and F^'^(X), respectively. This yields 



i-iLo + . k + i.)5f{p, K) = -ge.,V^F^-^{K)di^^f{p) + ^^(p) - f{p)) + ^^'-^ / Sf{p', K) , (Al) 



,V^F^^'^{K)d^J'^P{p) + *K/c^q(p) - /*(P)) + ^yfUv) (/p, 5r(p\ K)) /iVeq 

•J J yp, J — 

Defining 



SP (p, K) = —— — ' ' , ^"^ ^ ^ '- . ( A2) 



SfoiP^K) = [-ige,V^F^''iK)diP^rip) + t,.if^^ip) - f (p))) D-\K,yr,iy) , (A3) 
with D{K, V, i^) = Lu — Vi ■ v + iv we can write 

Sf\p,K)^SfUp,K)+z,yD-\K,^,,iy)l^ f 6f^ip',K) 

-''eq Jp' 



p" 



(A4) 



Using the shorthand notation 

r^{K)= I Sp^ip^K) (A5) 

and 

\{K, ^) - ^ / fl^{v)D-\K, V, v) (A6) 

-'"oq Jp 



we finally have 

<5f (p, K) = 5.ap, K) + ivD-\K, V, u)^-^^ri{K) (l + A + + . . .) 

which translates to the final result for the current by using 

Jp 

where we reintroduced the color indices. 
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